Large variation in the boundary-condition slippage for a rarefied gas flowing between 

two surfaces 
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We study the slippage of a gas along mobile rigid walls in the sphere-plane confined geometry 
and find that it varies considerably with pressure. The classical no-slip boundary condition valid 
at ambient pressure changes continuously to an almost perfect slip condition in a primary vacuum. 
Our study emphasizes the key role played by the mean free-path of the gas molecules on the inter- 
action between a confined fiuid and solid surfaces and further demonstrates that the macroscopic 
hydrodynamics approach can be used with confidence even in a primary vacuum environment where 
it is intuitively expected to fail. 

PACS numbers: 47.15.Rq, 47.10.ad, 47.61. Fg, 07.79.Lh 



It is traditionally assumed [l[ that in a fluid flowing 
along a solid surface, molecules nearest to the surface are 
globally stopped due to friction and collisions. This so 
called no-slip boundary condition has been very success- 
ful in modeling macroscopic experiments and it indeed 
forms one of the fundamental axiom of classical hydro- 
dynamics. However, it has recently been recognized that 
this standard condition is often not valid at submicro- 
and nanoscales 0, Furthermore, the hydrodynamic 
behavior close to a solid surface changes drastically with 
interfacial phenomena like roughness or surface chem- 
istry and the exact physics underlying these variations 
is not well understood [2-5]. Beyond its fundamental in- 
terest, elucidating these boundary conditions becomes a 
key issue for micro- and nano-electromechanical systems 
(MEMS-NEMS) such as sensors and actuators working in 
fluidic environments (liquid or gas). Although important 
results have been reported in a liquid environment (e.g. 
flow through nano- and microchannels or nano-tribology 
[3,SS,[3) there are only a few indications to the signif- 
icance of these phenomena in gases 

In this letter, we analyze a simple model apparatus 
able to continuously tune in the sphere-plane confined 
geometry the slippage boundary conditions at the solid- 
gas interface. By decreasing the surrounding pressure of 
a sphere facing a rigid wall in a gas (air or Helium), we 
find that these boundary conditions continuously evolve 
from a viscous regime supporting no slip to a ballistic 
regime with perfect slip. We interpret our results in terms 
of a giant modification of the gas slippage at the inter- 
faces. Therefore, our experiments appear to reconcile in 
a single setup boundary conditions that look confiicting 
at first sight. 

The experimental setup (Fig. [Ij is a homemade 
atomic force microscope (AFM) working at 300 K under 
controlled atmosphere. An optical-fiber based interfer- 
ometric detection of the cantilever [lit provides the 
required sensitivity to measure the impact of gas confine- 
ment on the viscous damping of the probe. The latter 
is a 460 /im long, 50 /im wide, and 2 /im thick, silicon 
AFM microlever with a i? = 20 /im radius polystyrene 
sphere glued at its extremity. In order to control the 
electrostatic interaction (see below), the whole probe 



radius (R) 



gap (z) ° 

gas pressure (P) 



FIG. 1: (a) Scheme of the experimental setup, (b) Zoom on 
the interaction zone between the facing surfaces defining the 
relevant physical parameters. 



(i.e. cantilever and sphere) is coated with a 200 nm 
thick gold layer, and the AFM chip is glued with sil- 
ver paint on a holder attached to the microscope frame. 
The probe spring constant k — 0.45 N/m and the reso- 
nance frequency /o = 9420 Hz have been measured using 
the Brqwnian motion of the thermally actuated lever at 
300 K 



12l |14| . The planar surface facing the sphere is 



a silicon substrate coated with a 200 nm thick gold layer 
and mounted on a high-precision positioning system to 
adjust the cavity gap z between the sphere and the sam- 
ple. An inertial motor makes sub-micron steps over a 
large 7 mm displacement range, whereas a piezo-scanner 
corrected for hysteresis distortions ensures a fine vertical 
positioning over a 1.5 /im range. 

For our measurements, it is important to conceive 
a method for measuring the absolute gap z. This is 
achieved by applying a voltage bias V to the probe with 
respect to the facing fiat surface as commonly done in 
Casimir force measurements where contact between fac- 
ing surfaces must be avoided [Hi U^. The cantilever 
is mechanically actuated with a dither at its resonant 
frequency using a phase-locked loop device and the res- 
onant frequency shift A/ = a {V — VcY is recorded as a 
function of V to extract the coefficient a that bears the 
desired distance information {Vc is the difference in ma- 
terials work functions). For z <^ R in the sphere-plane 
geometry, we use the relation a ~ /o7reoi?/(2fcz^) to find 
the absolute distance z. In fact, to reach a better preci- 
sion of 2% on z, we record a during a precise 1.5 scan 
aroimd the mean position and fit the results against the 
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FIG. 2: (a) Resonance spectra of the cantilever-sphere probe 
at three different pressures: ambient (Ibar), hehum gas 
(lOmbar), and vacuum (4 x lO^^mbar), when the fluid is 
not confined (the sphere to surface distance is z m 200 /^m). 
Note that the mechanical excitation has been reduced at low 
pressure to keep a similar amplitude at resonance, (b) Reso- 
nance spectra of the probe in vacuum for three different gaps 
z. The resonance width, i.e. damping factor, increases with 
decreasing separation. 
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FIG. 3: (a-c) Evolution of the damping factor T{z) = Fq -I- 
Th{z) as a function of the gap z between the sphere and the 
plane at three pressures: (a) ambient air at 1 bar, (b) helium 
gas at lOmbar, and (c) air vacuum at 4.10"^^ mbar. The sym- 
bols represent the data (at 1 bar and 10 mbar, the experimen- 
tal error bars are within the symbol size) and the continuous 
lines are the fitting curves according to the theoretical model, 
(d) Evolution of Th{z) with the gap z using the same color 
code as (a-c). The two curves (thin lines) surrounding each 
fitting curve (thick lines) provide an error estimate on the 
slip length b (from top to bottom: 6 = 0, 0.05, 0.1, 55, 70, 85, 
1000, 1200, 1400 Atm). 



above formula. 

We have measured the vibration amplitude A as func- 
tion of the excitation frequency / for three gas pressures 
P and several distances z. Examples of resonance curves 
are shown in Fig. [D^a) at large distance and for different 
pressures and in Fig. [^Jb) for different distances in a pri- 
mary vacuum. By fitting each curve with the Lorentzian 
response of an harmonic oscillator: 



A[f] - Ao 



7/0 



(1) 



we obtain the dissipation constant F — k^/ (27r/Q ) of the 
probe, with a fitting error smaller than 1%. The results 
are shown on the left graphs of Fig. [3] for three pressure 
conditions: (a) air at atmospheric pressure (1000 mbar), 
(b) helium at low pressure (lOmbar), and (c) air vacuum 
(4 X IQ-^mbar). 

When the sphere is far from the sample, the damping 
factor reaches a constant value Fq, characterizing the dis- 
sipation of the {cantilever -|- sphere} oscillating system in 
the fluid. This value is directly extracted from the data 
and reported in Table U (except in vacuum where the 
damping factor is not saturated at the largest distance 
and Fq was obtained by fitting with Eq.[3]as described be- 
low). As is also visible in Fig.[2lja), Fq decreases with de- 
creasing pressure P in agreement with previous works 

In the viscous regime this is mainly due to the exis- 
tence of a boundary layer of thickness (5_b oc 
representing a dissipation channel at finite frequency that 
adds to the intrinsic losses Fjnt of the lever (i.e. losses in 
the hmit /o ^ 0): Fo = Fi„t + 0{VP) Um\- 
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10 
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TABLE I: Measured slip length b and asymptotic damping 
rate Fq together with the calculated mean free-path Am and 
accommodation coefficient pd (see Eq. |4| at each pressure P. 



When the gap decreases, the hydrodynamic force 
due to the gas confinement contributes a z-dependent 
additional dissipation channel Th{z) such that F(z) — 
Fo + F/f(z). The central result of our work is shown 
in Fig. El^d) where Th{z) is plotted for different pres- 
sures. It is clearly seen that the dissipation observed at 
small distance is strongly reduced at low pressure. If in 
agreement with usual statistical mechanics [l[ we accept 
that the fluid viscosity rj does not depend on pressure 
(a reasonable assumption in the viscous regime) we can 
conclude that the boundary conditions at the solid-fluid 
interface should strongly change with pressure. In other 
words, the so-called slip-length b, usually used to charac- 
terize the fluid flow at the interface 0, 11, 0] , varies by a 
large amount, b is related to the fluid velocity gradient at 
the solid surface by surface = b surface (where v is 

the tangential fluid velocity) and can equivalently be in- 
terpreted as the fictitious depth below the surface where 
the no-slip boundary conditions would be satisfied. 

In the sphere-plane geometry with no-slip boundary 
conditions (i.e. 6^0), the dissipation constant Th{z) is 
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given in the hmit z <C -R by the Taylor formula [l|, Q : 



(2) 



where rj is the dynamic viscosity which does not change 
with pressure (1.8 x 10~^ kg.m~^.s~^ for air and 1.9 x 
10^^ kg.m^^.s^^ for Helium). To take into account the 
gas slippage at the boundaries, we follow the reasoning 
of Hocking [l^ and Vinogradova [19j and introduce a 
correction function /* such that 



(3) 



with J*{x) = 2x [(1 + x) In (l + |) - l] • In this for- 
mula (obtained for an incompressible fluid in the lami- 
nar regime po| ) we assume the same values for b on both 
surfaces since the two walls are coated with the same 
material. By fitting the three sets of data in Fig. |3Kd) 
with Eq. [3] in the range z < 20 /zm, we determine a slip 
length b for each pressure. Note that b is the only free 
parameter in the fit, except for the data in vacuum where 
both Fq and b are simultaneous adjustable parameters. 
b is however not sensitive to the choice of Fq since b is 
determined essentially by the absolute changes of F(z). 
The results are presented in Table H] with an estimation 
of the error shown in Fig. [31(d). The strong increase 
of b at lower pressure clearly shows that the friction of 
the confined fluid along the solid boundaries consider- 
ably changes with pressure going from the usual no-slip 
condition in ambient air (i.e. Th follows a 1/z law) 
to a quasi-perfect slip regime at low gas pressure (i.e. 
Th oc -ln(z/6)/&). 

In order to visualize the impact of a finite b on the 
fluid dynamics close to the surface, we compute the (ra- 
dial) fluid velocity proflle v{Z) in the gap z between the 
two surfaces using the analytical theory of 3 Due 
to the fluid incompressibility and the limit z <^ R, the 
fluid is essentially ejected from the gap in the radial di- 
rection in response to the vertical displacement of the 
sphere 2l|. The comparison between the three differ- 



ent values of b obtained in the experiment (see Table |T| 
is shown in Fig. SI We clearly switch from the usual 
Poiseuille parabolic velocity proflle (i.e. ?; ~ on the 
solid boundary) for 6 — > 0, to a quasi-constant velocity 
profile in the gap for b — )■ +oo. Therefore our experi- 
ment reveals a continuous transition between these ex- 
treme regimes. 

Now, from a microscopic point of view, the slippage 
coefficient is linked to the very nature of the interaction 
between moving surfaces and air molecules. Following 
the historical approach of Maxwell's kinetic theory of 
gases j22!j, two interaction channels can be distinguished: 
a specular one where molecules collide elastically with the 
surface and a diffusive one where molecules are reflected 
diffusively by the wall. The slip length b in this statistical 
model is given by the Maxwell formula IJ, |22|424| : 
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FIG. 4: The normalized radial fluid velocity profile v{Z, r)/vo 
(«o is the sphere velocity) as function of the coordinate Z in 
a 2: = 1 fim gap and for a radial coordinate r = 1 /^m. The 
profiles are calculated for the three b values of Table IT] 



where Xm is the typical mean free-path of gas molecules 
and Pd the tangential momentum accommodation coeffi- 
cient, i.e. the fraction of those particles hitting the sur- 
face with a diffusive reflection. 

The main dependence of b with pressure comes from 
the molecular mean free-path (see Table [l|: 



3 Pd 



(4) 
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where T is the temperature, ks the Boltzmann's con- 
stant, and a the molecular cross-section. Two important 
asymptotic regimes are clearly found in the experiment 
(see Fig. d]). On the one hand, when the mean free- 
path is extremely small compared to other macroscopic 
dimensions (i.e. at ambient air pressure), the fluid par- 
ticles interact strongly with themselves, and even more 
strongly with surfaces located within a characteristic dis- 
tance given by Am- This is a consequence of multiple 
collisions, reflections, and adsorption processes 0-0, [l^]. 
In this diffusive regime, the tangential velocity of the 
molecules decreases at short distance from the surface, 
such that the slip-length b tends to zero and the no- 
slip condition applies. On the other hand, when the 
mean free-path is considerably larger than the gap, the 
molecules interact mainly with surfaces and there is no 
momentum transfer among the molecules themselves. In 
this ballistic or molecular flow regime, the velocity gra- 
dient vanishes at the surface. This results in a giant slip 
length b. 

In Eq.[4l the accommodation coefficient pd (see Table|I]) 
depends on the surface properties and gas density (i.e. 
on the probability of multiple collisions). Recent analy- 
sis based on the fluctuation-dissipation theorem and the 
Green-Kubo relation emphasize the importance of sev- 
eral microscopic parameters such as surface roughness 
and surface defects on the molecular dynamics close to 
the surface In our experiment, the rms roughness of 

the gold surfaces was found by AFM to be i5r ~ 3 nm for 
both the sphere and sample. This certainly contributes 
to increase pd and reduce b as compared to a perfect sur- 
face. 

Note that Eq.[3]is an analytical solution of the Navier- 
Stokes equations whose validity is well established when 
the fluid density is sufficiently large to provide a local 
equilibrium within the fluid. However, in the present 
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experiment in vacuum, the mean free-path Am is much 
larger than any other relevant length (e.g. the gap z) so 
that the system enters in a molecular flow regime which 
should be discussed within the more general frame of 
Boltzmann's kinetic equations [H, [H- Therefore, the slip 
length b deduced from our data at low pressure should 
be considered as an extrapolation showing the limita- 
tion of the usual hydrodynamic approach. Because the 
low-pressure data of Fig. [3] are well reproduced in our 
analysis, our work stresses the remarkable robustness of 
the Navier- Stokes Eq. |3]even in this molecular regime. 

Finally, it is worth commenting on the differences be- 
tween the present sphere-plane experiment and recent 
results obtained in a plane-plane geometry [Ij]. At at- 
mospheric pressure, we obtain here a small slip length in 
agreement with Refs. d, [3| whereas Ref. [15] reports on 
a perfect slip in the same fluid (air) . One of the difference 
is the geometry, which is indeed known to be a critical 
parameter at micro- and nanoscales BSEHIii. An- 
other difference in Ref. is the probe velocity which 
resulted from the thermal motion of the cantilever with 
oscillation amplitude A 0.05 nm, whereas the probe is 
here mechanically actuated i.e. A « 30 nm (in the same 
context compare Refs. [l3| and [1]). In addition, the sur- 
face roughness was only a few Angstroms, resulting in a 
smaller accommodation coefficient pd and a larger slip- 
length b. 

In conclusion, we have discovered that the slippage of a 
gas along mobile rigid walls varies considerably with pres- 
sure in the sphere-plane confined geometry. The classi- 
cal no-slip boundary condition valid at ambient pressure 
changes continuously to an almost perfect slip condition 
in vacuum. This study emphasizes the key role played 
by the mean free-path Am on the interaction between a 
fluid and solid surfaces and demonstrates that the macro- 
scopic hydrodynamics approach can be used with confi- 
dence even in good vacuum conditions. We anticipate 
that our work will have an impact on the MEMS and 
NEMS engineering and will motivate further fundamen- 
tal studies of the physics of gas slippage along solid and 
mobile surfaces. 
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